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(^ We present holographic descriptions of thermahzation in confornial field theories 



using probe D-branes in AdS x S space-times. We find that the induced metrics 
on Dp-brane worldvolumes which are rotating in an internal sphere direction have 
horizons with characteristic Hawking temperatures even if there is no black hole in 



^ the bulk AdS. The AdS/CFT correspondence applied to such systems indeed reveals 

thermal properties such as Brownian motions and AC conductivities in the dual 
conformal field theories. We also use this framework to holographically analyze time- 
dependent systems undergoing a quantum quench, where parameters in quantum 
field theories, such as a mass or a coupling constant, are suddenly changed. We 
confirm that this leads to thermal behavior by demonstrating the formation of 
apparent horizons in the induced metric after a certain time. 
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1 Introduction and Summary 



The AdS/CFT correspondence [l]-^ has been useful in understanding time-dependent 
processes on both sides of the duahty. On the gravity side, some insight has been obtained 
for time-dependent gravitational backgrounds like space-like singularities using the gauge 
theory dual |5]-[7]. However, the correspondence has been most successful in the explo- 
ration of non-equilibrium properties of strongly coupled field theories at finite temperature 
and/or finite density using their gravity duals. The significant results in this area include 
the well known predictions for linear response functions [8 11 . In recent years, however, 
this approach has been extended to the non-linear domain as well, e.g. in the derivation 



of non-linear fluid dynamics from gravity 12 and in the study of thermalization in terms 



of the dual process of black hole formation 13 14 



A particularly interesting class of non-equilibrium problems involves time-dependent 
couplings or masses which pass through a critical point either at finite temperature or 



at zero temperature 15 . This problem has experimental relevance in cold atom physics 



and has been studied theoretically in two extreme limits. In the first limit, a parameter 
is changed quickly through a critical point (a quantum quench) and the subsequent time 
evolution is studied 16 22 . In the other limit, the parameter changes slowly across a 



critical point 23 . In both cases, dynamical quantities have universal signatures of the 



critical point. There are, however, few conventional analytical tools to investigate such 
problems. It would be clearly interesting if gauge-gravity duality can be applied to such 
situations. 

In this paper we take a small step in this direction. We investigate the dual bulk 
dynamics of a class of conformal field theories with time-dependent couplings, starting 
in an initial vacuum state. In known field theories, this is an interesting problem even 
when the couplings do not pass through a critical point. For example, if we start with a 
free massive field theory in its vacuum state and suddenly change the mass to a different 



value, the resulting real time correlation functions become thermal at late time 17,19 



Similar behavior can be established in some interacting field theories as well 19,22 , but 



the tools for investigating such phenomena in strongly coupled theories are rather limited. 
When the field theory has a gravity (as opposed to full fledged string) dual, a time- 
dependent coupling corresponds to exciting a non-normalizahle mode of the dual bulk field 
and calculating the time evolution. This is quite similar to the investigation of cosmologi- 



cal backgrounds in 5 , forced fluid dynamics 24 or black hole formation 13,14 . However 



the class of field theories we consider are considerably simpler, and may have different 
applications. These are defect or flavored conformal field theories resulting from probe 



D-branes in a background space-time of the form of AdSm x "S*". Such probe branes 25 



have been used to model flavor physics 26-28 and quantum critical phenomena : this 



"top-down" approach 29 33 is complementary to the "bottom-up" approach involving 
the physics of charged black holes in AdS [lO||lT|[34|[35]. One positive feature of this 
top-down approach is that in this case the dynamics of the probe branes describe the 
strongly coupled dynamics of a known field theory. Furthermore, since the gravitational 
background is not altered in the probe limit, the entire dynamics is given by the world- 
volume dynamics of the probe brane. Solving this dynamics is usually easier than solving 
Einstein equations. 

In the following, we will find nontrivial time-dependent classical solutions of several 
kinds of probe branes which represent time-dependent couplings of the dual defect confor- 
mal field theories. We find that the induced worldvolume metric on a large class of such 
classical solutions have apparent horizons characterized by a temperature in spite of the 
absence of black hole horizons in the AdS bulk metricQ In fact, small fluctuations of the 
worldvolume fields around these solutions behave like fields in space-times defined by this 
induced metric. From the point of view of the dual field theory, this means that an initial 
vacuum state can evolve into a thermal state in the presence of time-dependent couplings. 
The emergence of apparent and event horizons on the worldvolume is somewhat similar 
to the emergence of an acoustic horizon for phonons around a supersonic fluid flow 1431 . 
Indeed, this kind of phenomenon can occur in a large class of non-linear field theories [4^ . 

It is very important that in our probe approximation the dual field theory should 
be regarded as an open system. Consider such a setup in AdS^ x S^ dual to the four- 
dimensional TV = 4 super Yang-Mills (SYM) which couples to a defect or flavor conformal 
field theory (see Fig. [I|. Though the energy supplied by the time-dependent coupling first 
excites the defect or flavor sector, this energy will finally dissipate into the A/" = 4 SYM 
sector due to the interactions between each other. In the bulk description this means 
that energy can flow from the probe brane to the bulk gravitational degrees of freedom. 
Indeed we find that typically the worldvolume apparent horizon evolves into an event 
horizon similar to that of a static black hole space-time only when the time derivative of 
the classical solution approaches a constant. In such a case, the long time correlators of 
fluctuations become thermal with the Hawking temperature of the event horizon. This 



■' Here we would like to mention earlier works on the presence of horizons in the induced metric. A 
general analysis of appearance such horizons on accelerated D-branes has been done in [36| . World-sheet 



horizons of F-string has already been discussed in 37-40 . In 41 , the presence of horizons on time- 



dependent D-branes in flat space-time has been analyzed. In 42 horizons on D-branes from accelerated 
observers have been studied in pure AdS space-time. 



can be explicitly demonstrated for the cases where the fluctuations can be decoupled. For 
most cases, these fluctuations are coupled in a non-minimal way. However, one expects 
that correlators will still inherit the thermal behavior implied by the background metric. 
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Figure 1: A schematic diagram for our non-equilibrium systems in the AdS/CFT. Here 
we concentrate on the example of a probe D7-brane in AdS^ x S^ discussed in section 4. 
The other cases can be interpreted in the same way after obvious modiflcations. 

In summary, our holographic setups can be regarded as an example of non- equilibrium 
steady states. One of the most important non-equilibrium properties of our systems is that 
there are two different temperatures: one is the bulk temperature Te which is vanishing 
for pure AdS and the other one is the Hawking temperature T^ of the apparent horizon 
of the induced metric of probe D-branes. 

The simplest example in our setups is a probe D-string or F-string extended along the 
radial direction of AdS^. In the dual A/" = 4 SYM theory this represents a monopole or a 
quark. In this case the relevant class of classical solutions can be found analytically. In 
these solutions, the end-point spins along a direction on the S^ leading to a wave moving 
along the string towards the Poincare horizon. In the dual defect CFT this corresponds to 
a certain time-dependent coupling of the hypermultiplet fields with M = 4 vector multiplet 



fields 45,46 . This generates a R-charge chemical potential only for the hypermultiplet 
and not for the N' = 4 SYM sector. If we start with a static string in the far past and 
turn on the spin, an apparent horizon develops on the worldsheet. When the value of 
the spin approaches a constant at late times, the induced metric coincides with that of 
the BTZ metric with an angular coordinate suppressed. The Hawking temperature of the 
event horizon is proportional to the final angular velocity u, given by simply Th = ^■ 
Fluctuations around this classical solution satisfy simple equations in this background 
metric leading to correlators which are thermal. 



In particular, this leads to a Brownian motion of the end point. In the dual field 
theory this manifests itself as a Brownian motion of the monopole position as well as a 
Brownian motion in the internal space 5*0(6). When the string extends all the way to 
the boundary, the mass of the monopole or the quark is infinite. However terminating 
the string on a suitable space-filling D7-brane makes the mass finite. It should be again 
emphasized that this is a non- equilibrium situation : it is only the monopole or quark 
which behaves thermally while the other A/" = 4 degrees of freedom are still at zero 
temperature. 1/A^ corrections will eventually thermalize the latter - but that requires 
going beyond the probe approximation. Since we have an open system, energy can fiow 
into the vector multiplet sector. In a closed system, we typically expect thermal behavior 
if the time dependence vanishes at late times. In contrast, in our case we need a steady 
time dependence to achieve thermality since energy can fiow out of the defect CFT at a 
constant rate. Such Brownian motion behavior has been observed for fiuctuations around 
static strings in the background of AdS black holes, i.e. when the field theory is already at 



some finite temperature 47,48 . In our case, the parent field theory is at zero temperature 
- an effective temperature is produced by a time-dependent coupling. 

The physics of other probe D-branes is quite similar, but requires a combination of 
analytical and numerical tools. For D3, D5 and D7-branes in the AdS^ x 5*^ geometry, 
we find classical solutions with angular momenta dual to R-charges which lead to induced 
worldvolume metrics with horizons. For a probe D3-brane, which is dual to a (2 + 1)- 
dimensional defect CFT, we calculate the conductivity by turning on a worldvolume gauge 
field. There is a finite AC conductivity which arises from dissipation into the N' = 4 de- 
grees of freedom. A similar behavior has been observed earlier using probe D-branes at 



finite temperature, i.e. in AdS black hole backgrounds 29,33 . In contrast, we observe 
this behavior in a zero temperature field theory where a time-dependent coupling effec- 
tively produces a thermal state. In the D7-brane case, we can rotate a D7-brane in the 
standard D3-D7 system [26] dual to four-dimensional A/" = 4 SYM with a massive fiavor 



hypermultiplet. This setup has been already discussed in |!49| (see also 50 for a similar 
analysis for D3-D5 system). This rotation excites a R-charge only in the hypermultiplet 
sector. We find that for a sufficiently large R-charge, the induced metric of the probe 
D7-brane has a horizon which surrounds a cusp-like singularity on the worldvolume. This 
structure looks very similar to the structure of the standard black holes in general rel- 
ativity. As a final example, we holographically describe a time-dependent mass in the 
defect CFT sector dual to a probe D5-brane. We numerically show that an apparent 
horizon is formed under a time-dependent change of mass, thus describing thermalization 
of the dual CFT. In free field theory, it is known that the effective temperature after the 



quench depends on the momenta of each modes since they are decoupled with each other 



in free field theory 19,21 . However, our holographic result predicts that this picture will 
be completely changed for the strongly interacting gauge theories. In the latter case the 
different momentum modes observe a common temperature. 

It is also an intriguing problem to associate an entropy to such a D-brane with thermal 
horizon. We will show that one holographic candidate for such a entropy is divergent in 
general except in the lowest-dimensional case i.e. Dl-brane case. We suggest that this 
might be interpreted as an entanglement entropy when we trace out the Af = 4 SYM 
sector, though we would like to leave its detail for future work. 

This paper is organized as follows: In section 2, we consider probe Dp-branes which 
are extended in AdS space-time and rotating in a sphere direction. We show that an 
apparent horizon is always formed in the Poincare AdS space. We calculate the energy 
dissipation in such systems. In section 3, we study the thermal property of the rotating 
Dp-brane solutions constructed in section 2. In particular, we will study Brownian motion 
and thermal correlation functions for a probe Dl-brane and the AC conductivity for a 
D3-brane. In section 4, we analyze the rotating D7-brane with a mass for flavor hyper- 
multiplet. We will show that for a sufficiently fast rotation, the induced metric of the 
probe D7-brane has a horizon. In section 5, we will give a brief analysis of a holographic 
dual of a quantum quench caused by a rapid change of a mass parameter in a defect CFT. 
We consider a particular case of probe D5-brane and numerically show the appearance of 
an apparent horizon in the time-dependent induced worldvolume metric. 

2 Worldvolume Horizons and Mini Black Holes from 
Rotating D-branes in AdS 

Typical setups of AdS/CFT correspondence are AdSd+2 x -^^ backgrounds in string 
theory or M-theory, where AdSd+2 is the {d + 2)-dimensional AdS space in Poincare 
coordinates and X'^ is a g-dimensional Einstein manifold fllfl]- In this paper we mainly 
assume the simplest case X'^ = S"^. We can write their metric as follows (we set the AdS 
radius R = 1 throughout this paper, for simplicity of presentation.) 

ds^ = -r^de + r^Y^ dx] + ^ + {d9^ + sin^ 9d^^ + cos' 9dnl_^) . (2.1) 

This clearly includes the celebrated example of the AdS^ x S^ in type IIB string theory 
by setting d = 3 and q = 5, which is dual to the four-dimensional J\f = 4 super Yang- Mills 
theory. 



The D-brane systems we would like to consider now are Dp-branes rotating in the 
S'^ directions. These rotations are dual to non- vanishing R-charges in CFT^+i. The 
equations of motion require that such D-branes should rotate along the largest circles in 
S"^. Below we will show that such D-brane systems can be regarded as mini black holes 
in AdS spaces and they describe locally thermal vacua in a probe limit. 

Since the analysis of rotating Dl-branes is much simpler we will first study this and 
then proceed to the analysis of general Dp-branes. 

All the rotating solutions that will be constructed in this section and their relation to 
dual gauge theory are summarized in Fig. [2j 
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Figure 2: We present the schematic pictures of (a) a thermal Dl-brane solution extended 
in Poincare AdS^, (b) a suspended Dl-brane solution in Poincare AdS^, (c) a Dl-brane 
solution in global AdS^, and their CFT duals: (a'),(b') and (c'). The Dl-branes are all 
rotating along a great circle in the S^. 



2.1 Rotating Dl-branes 

Consider a probe Dl-brane (or a F-string equivalently owing to S-duality) in AdSd+2 x S'^ 



(2.1). The Dl-brane extends in both t and r direction, infinitely and locahzed at Xi = 0. 
In addition, it is spinning in the ip direction. Therefore, the Dl-brane world-volume is 
specified by the function 

ip = ip{t,r), (2.2) 

at constant 6{= 9^). 

In this setup, the DBI action of the Dl-brane is 

Sdbi = -Tdi / dtdrL , (2.3) 



sin^/9 

L = \l 1 + r^^'y sin^ 9, ^0^ , (2.4) 



where Tdi = ^^^, is the tension of a D-string. The equation of motion reads 

d f r"^ siv? 9^'{)'\ d /'sin^6'^,0' 



dr \ L J dt \ r'^L 
Consider solutions of the form 



(2.5) 



/I -uj'^lr'^ 
^ = ujt + g{r), ^'(r) = ^-^-^, (2.6) 

0. = l. (2.7) 

where A is an arbitrary constant. The Lagrangian density for this solution is 



In order to avoid a negative sign inside the square root, we need to require 

^=1. (2.9) 

In this way, we find the rotating Dl-brane is given by the solution 

^ = iot-- + ^Q, (2.10) 

r 

up to an integration constant y^o- 

The induced metric on this Dl-brane worldsheet is 

ds^ = -{r^-UJ^)de + 2^dtdr+(-^ + '^\dr\ (2.11) 



We can rewrite this by defining a new time coordinate 

r.t-1-^logl^. (2.12) 

as follows 

^^2 ^ _^^2 _ ^2)^^2 ^ ^!:1^ _ (2.13) 

Notice that r becomes identical to t as we approach the boundary r — )■ oo. 



Interestingly, the induced metric (2.13) coincides with the BTZ black hole with the 
angular coordinate suppressed. Therefore it has a horizon at r = cu and we can read off 
its Hawking temperature Tu 

Th = ^, (2.14) 

ATI 



by Wick-rotating r into a Euclidean time (see the picture (a) in Fig. |2]). 

If we take into account the back-reaction of this solution to the supergravity back- 
ground, it is natural to expect that such Dl-branes produce a very small black hole in 
the bulk AdS localized in the R^ direction (we may call it a mini black hole). 

This shows that the rotating Dl-brane describes a thermal object with temperature Th 
in the dual CFT . Notice that the bulk of AdS is at zero temperature. If we concentrate on 
the most important example of AdS^ x S^, our system is dual to A/" = 4 super Yang-Mills 
theory coupled to an infinitely heavy monopole. The A/" = 4 gauge theory itself is at zero 
temperature, while the monopole is at finite temperature T^. Therefore such systems are 
in non-equilibrium steady states. 

The rotating Dl brane solution corresponds to a time dependent coupling in the 
A/" = 4 theory coupled to hypermultiplets living on the zero dimensional defect. The Dl- 
D3 system is 1/4-BPS and the D1-D3 open strings lead to the two complex scalars (Q, Q) 
of hypermultiplets which belong to the fundamental and anti-fundamental representations 
of the color SU{N) gauge group. Let us express the three complex adjoint scalar fields in 
the Af = 4 super Yang-Mills by ($i, $27 ^3)- These correspond to cartesian coordinates 
in the transverse C^ composed of (r, ^5) where Q5 represents the 5-sphere. We choose 
$3 such that its phase rotation describes the one in the (p direction and that 6 = -r 12 is 
equivalent to $1 = $2 = 0. We now argue that the time dependent coupling term which 
corresponds to a uniformly rotation Dl-brane is given by 

I dt\Q [Im(<l>3e-*"*)]' Q^Q [lm($3e-*'^*)]' q] . (2.15) 

To see this consider going to the Coulomb branch where ($3) 7^ 0. This corresponds 
to separating out some number k of the original A^ D3-branes. Apart from generating 

9 



masses for the off diagonal vector multiplet fields, the coupling in (2.15) generates masses 



for appropriate components of the hypermultiplet fields. Consider the simplest situation 
where only one D3-brane is separated in the 6 = 7r/2 hyperplane, so that we only have 
($3) = diag[$3, 0, 0, ■ ■ ■ ]. In this case the topmost component of Q would become massive. 
This mass should be equal to the length of the shortest open string joining the Dl-brane 
with this separated brane (upto a factor of the string tension). The rotating Dl described 



iip{t) 



ze 



iuit^ 



in the transverse C^, where z takes any real 



above is described by (0, 0, ze 
value. The shortest distance is then given by |Im($3e~*'^*)|, leading to a mass of the 
scalar fields [Q, Q) exactly as predicted by (2.15). In the above discussion we have taken 
($1) = ($2) = to extract the main time-dependent part of the interaction. However the 
result can be generalized to a generic point in the Coulomb branch, where other terms 
involving the $1 and $2 would be involved. 



As a further check on consistency, note that the operator in (2.15) has dimension 1. 



As we will see in a later section the field (f is a massless minimally coupled scalar on the 
worldsheet. The usual AdS/CFT formula (now applied to AdS2 on which the Dl-brane 
wraps) relating bulk masses to conformal dimensions then predict that the dual operator 
should have dimension 1. 

We thus see that the AdS/CFT correspondence shows that if we consider a system 
with such time-dependent interactions, the system gets thermal (see Fig. [I]). 

2.2 General Fl/Dl solutions 

It turns out that there are more general purely right-moving or purely left-moving solu- 
tions of Dl or Fl branes in the background of the class of metrics which include AdS x 5* 
as well as black branes or black holes in Poincare or global coordinates. Consider a metric 
of the form 



ds' 



-f{r)dt' + 



dr"^ 

W) 



Y^ dxl + {de'^ + sin^ Qdip^ + cos^ ^rffij. 



2) ) 



(2.16) 



j=i 



and introduce the outgoing and ingoing Eddington-Finkelstein coordinates 

dr r dr 



fir) ' 



u = t ■ 
The DBI action then becomes 

S = —Tdi I dudv f{r) 



v = t + 



m 



(2.17) 



fir 



-du^dyif 



1/2 



(2.18) 
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The equation of motion then becomes 

dudyif + -dyifidu I . . 1 + -Ou^Oy I , . ) = ' (2-19) 

where 

L = l- -—duipdyip . (2.20) 

fir) 



It is easy to see that any Lp which satisfies either du^ = or d^ip = is a solution of (2.19). 
In particular, solutions which are functions of v alone may be thought of the retarded 
effect of a boundary value of (p. 

The induced metric produced by such a retarded solution is given by 

(is2^^ = -f{r)dudv + {dy^fdv"^ = 2drdv - [/(r) - {dyipy]dv^ . (2.21) 

This is a two-dimensional AdS Vaidya metric which shows the formation of an apparent 
horizon at /(r) = [dyipY, provided this equation has a solution for real r. The dual 



CFT description is captured by the interaction (2.15) with ut is replaced with the general 



function (p(t) as is clear from the derivation in (2.15). 



Special cases of the background metric (2.16) include Poincare patch AdS^, with 
/(r) = r^ and black 3-branes with /(r) = r^(l — ^). With the replacement of Ylii^^iY ~^ 
dfll this can include global AdS^ as well as AdS-Schwarzschild. We will discuss these so- 
lutions further in the following sections. 

In the special case where the function p>{v) approaches p>{v) — > uv at late retarded 
times, the apparent horizon asymptotes to an event horizon at /(r) = w^ which signals 
thermalization. In dual CFT, this system consists of two sectors with different tempera- 
tures. 

2.3 Rotating Dj9-branes 

We can extend the previous system to the rotating higher-dimensional Dp-branes. Con- 
sider a rotating Dp-brane for p < d + 1 in AdSd+2 x •5''' space-time which extends in 
{t,r,xi . . . ,Xp-i) directions of AdS and rotates in ip direction at ^ = 7r/2. We again 
assume ip depends on only t and r directions. 
The DBI action is given by 



Sdbi = -Tp f dtdrd^'^x rP-^Jl + (rV' - ^) • (2.22) 



11 



It follows that we obtain the equation of motion: 

8, { , '-^"^ \ - 9. ( , ^"•^' I . (2.23) 

We now consider the rotating solution with angular velocity u along ip direction, and we 
thus assume the following form of the solution: 

ipit,r)=Lot + gir) . (2.24) 

Substituting the ansatz into the equation of motion the function g{r) looks like 



where A is the integration constant, and should be set to A = u^ to avoid a divergence 
of the function g{r) aX r = u. 

The induced metric on the rotating Dp-brane is 

ds'^i^^ = -(r^ - uj'^)dt'^ + (^ + 9^)dr'^ + ^ug'dtdr + r'^dxl_^ 

= -{r'- u')dT' + ^^^^/r' + rHxl_, , (2.26) 



where the new time coordinate r is defined as 



w^+i 



T = t- I dr ^^^^^^^^ . (2.27) 



This metric (2.26) clearly has a horizon at r = w. After Wick rotation of r we can 



determine its Hawking temperature by demanding smoothness at r = w: 

T = ^^— . (2.28) 

Itx 

2.4 Back-reactions and Charge/Energy Dissipations 

In the above analysis, back-reaction to the supergravity sector has been neglected since 
we employed the probe approximation. We would now like to see how much and when 
this is justified. Since the essential issues are the same, we concentrate on a rotating 
Dl-brane located at Xi = 0. 

The total energy of our rotating Dl-brane conjugate to time t is given by 

E= / dr ^ / dril + — ). (2.29) 

Jo /i + ^2(^/)2__^02 Jo \ r J 

12 



This shows that the energy density blows up at r = 0. Thus in the IR hmit we cannot 
neglect the gravitational back-reaction. This is intuitively understood if we plot the form 



of Dl-brane profile (2.10) at some given time in the (rsm(p,rcos(p) plane, as in Fig. 



|3| Notice that the Dl-brane wraps infinitely many times at r = 0. Similar behavior is 
present for the rotating D-brane in AdS black holes. 




Figure 3: The profile of a rotation Dl-brane in AdS at some given time in the (x^ = 
rsin(p,x^ = rcos(/9) plane. The time evolution is obtained by rotating the profile. 

It is reasonable to conclude that the large back-reaction will result in the formation 
of a black hole in the bulk of AdS, centered at Xj = r = 0. The size of this black hole 
should grow as the energy is pumped into it from the Dl-brane constantly. To obtain this 



energy flux we need to calculate the energy-stress tensor 51 for the probe Dl-brane. We 
define the energy-stress tensor by 



6S 



-9^9c 



-9ch 



(2.30) 



where gab denotes the bulk metric (the index a and h runs all ten (or eleven) coordinates 
of our space-time). This satisfies the equation of motion 



VaT\ = 



(2.31) 



which is reduced to 



daiT'^t^c 



(2.32) 
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for static space-times. This leads to the conserved energy 

E= I dry/^T\ . (2.33) 



In our rotating Dl, the exphcit expressions for T^ are 

, ^ = Tbi 1 + - 



t _rj. 1 + rV' _^ (^ ^ UJ' 



Thus we can calculate the time evolution of the total energy as 

Jdr^^T\ = JdrdriV^Tl) = ,/^T[ IZ"^ = Tmu'-To^u' = . (2.35) 



dE _ d 
dt dt 



Even though the total energy does not depend on time, this result shows that the energy 
'^^ — Tdiuj"^ per unit time is injected at the boundary r = oo by some external system 



dt 



reservoir) and the same energy -rr = TdiLO is dissipated from the IR region into the 



bulk AdS as summarized in the picture (a') in Fig. M This dissipation from the Dl-brane 
to the bulk AdS will create a localized black hole in AdS. The dissipation via the AdS 
horizon is very similar to the phenomenon called drag force, which is dual to the jet 



quenching 52 54 



One may wonder how we can explicitly understand the external injection of energy in 
our stationary solutions as we argued. To see this, let us introduce UV and IR cut off 
into our rotating Dl-brane system so that it extends from r = rjji to r = ruvi^ fiR)- 
In this setup it is clear that the energy flux Toioo"^ is coming from r = rjjv and the same 
energy is going out at r = r/^. The presence of this energy flux is clear from the non- 



vanishing value of T[ in (2.34) at r = r/^j and r = rjyy. At r = r^/j, the energy is not 
reflected back but its back-reaction will make a localized black hole, which absorbs the 
injected energy. Another way to see the energy flow explicitly is to consider general time- 
dependent solutions presented in section 2.2. In this case we can choose the Dl-brane 
profile such that the energy is injected like a pulse and this travels into the IR region. In 
this case, T[ gets non- vanishing at the UV cut off only during a certain time interval. Our 
stationary rotating Dl-brane solution can be regarded as an infinitely many collections of 
such pluses and therefore it requires the constant energy injection at the AdS boundary. 
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Similarly we can calculate the dissipation of angular momentum which is dual to a 
R-charge in the CFT. The angular momentum of a rotating Dl is given by 

-— = Tdi / dr— 
dip J r^ 



Qr = —.=Tdi I dr- , (2.36) 



which is divergent like its energy. 

The emission toward the horizon r = for the Poincare AdS can be estimated as 



dQ 



R 



dt 



-Tdi ^ 



r=0 



1 + r'^ip''^ — 02/r^ 



r=0 



-Tdico . (2.37) 



This amount of angular momentum (=R-charge) comes from the boundary r = oo and 
flows into the horizon r = 0. For higher-dimensional rotating Dp-branes, the situations 
are almost the same. From the viewpoint of the boundary CFT which lives on R^''^, only 
the codimension d — p + 1 part, defined by the intersection between the AdS boundary 
and the Dp-brane, is thermalized at the temperature Th, while the other part is at zero 



temperature in the probe approximation. The energy dissipation we found in (2.35) 
corresponds to the energy flow and the R-charge flow from the thermal part to the cold 
region in the bulk of R^''^. Therefore the temperature of the bulk region should constantly 
increase. This is dual to a growing black hole once back-reaction is taken into account. 

2.5 Rotating D-branes in Global AdS 

One way to suppress a large IR back-reaction is to replace the Poincare AdSd+2 with the 
global AdSd+2- Again we will concentrate on Dl-branes just for simplicity. Then we need 



to take /(r) = r + 1 in the equation (2.16). The rotating Dl-brane is described by the 
proflle 

(fi = uj{t + arctanr) . (2.38) 

Notice that the Dl-brane in the global AdS pierces the center of the AdS and thus has two 
end points at the boundary, which are identifled with the north-pole and south-pole of 
the boundary S^, respectively, as shown in the picture (c) in Fig. pi This is conveniently 



described by extending the range of r in (2.38) to — oo < r < oo. 



The induced metric on the worldsheet is given by 

dsld = -{r' + 1 - u;')dT' + ^, ^^['_ ^, , (2.39) 

where t = t + h{r) and h'{r) = (i,r'2)(%+i-u]^) - Therefore, a horizon exists only when 
|a;| > 1. This gap is analogous to the Hawking-Page transition between the global AdS 
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and the AdS black hole 55 . In our case, it is dual to the confinement /deconfinement 



phase transition of the hypermultiplets from D3-D1 strings. The Hawking temperature is 
given by Th = ^^^J? • 

The total energy of the Dl-brane is given by 

Its divergent part which is independent of u can be canceled by the standard holographic 
renormalization. The total angular momentum becomes finite for the Dl in global AdS 
as given by 

dr^— = ttTdiUJ . (2.41) 

Thus we can conclude that the back-reaction of the probe Dl-brane to the metric is not 
large even when u is non-zero. This is rather different from the rotating Dl-brane in 
Poincare AdS. 

Just as in the Poincare patch solutions in the previous sections, there are general 
purely left moving or purely right moving solutions. The solutions corresponding to the 
retarded effect of specified boundary conditions are of the form 

y} = y9(t + arctanr) , (2.42) 

for arbitrary choice of the function (f{v). The R-charge current shows the behavior 

^ = {r' + l)^Too = ^{t + vr/2) - 0(t - n/2) ^ j, - 32 ■ (2.43) 

This means that in the dual CFT, the current ji starts from the north pole Pi, it fiows 
into the bulk of 5*^ and is eventually absorbed at the south pole P2 as the current J2 
(see the picture (c') in Fig. [2]). The time delay At = tt is just the propagation of this 
current at the speed of light from Pi to P2. Notice that in this global AdS, there is no 
energy dissipation toward the bulk (i.e. TV = 4 SYM sector) at planar order, which will 
be because the bulk system is in the confinement phase. The current is carried by fiowing 
R-charged particles. Also notice that the effective temperature for P and Q are |^ and 

27r' 

2.6 Suspended Rotating Dl-branes in Poincare AdS 

As we have seen, in the global AdS space, the rotating Dl-brane always has two end 
points at the boundary. This is because we cannot have a single source of charge current 
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in a compact space. In the Poincare AdS case, this is no longer true since the space is 
non-compact. However, still it is interesting to consider a system where there are two 
point sources Pi and P2 in the boundary R^'^. In this case, there are two possibilities 
of rotating Dl-branes. The first one is two disconnected rotating Dl-branes infinitely 
stretching in the radial direction. As we already analyzed, this has a thermal horizon at 
r = u and has the temperature T = ^. We can take the temperatures of two points 
independently, denoted by Ti and T2. Then the total R-charge current flow from Pi to 
P2 in the dual CFT is given by 

dQn 



dt 
and the energy flow is 



27rTz5i(Ti-T2) , (2.44) 



^ = 47r^T^i(Ti^-T|). (2.45) 

At the two points, the system gets thermalized and the carrier possesses the energy u per 
a unit R-charge. They will strongly interact with the bulk gauge theory even though it is 
at zero temperature. Notice that in this case there is no correlation between two points 
and the physics does not depend on the distance L between them. 



We can calculate the R-charge conductance from (2.44). We consider the setup of 
type IIB string on AdS^ x 5*^ and recover the radius of the AdS space R = {AngsNa'^) ^ . 
The chemical potential at the two points are identifled with ui and a;2. Thus we flnd the 
conductance, which is the ratio of voltage V to the current /, is given by 

' «« To.R'^^^, (2.46) 



V Ui -UJ2 VTT a/a 

where A = Nqs is the 't Hooft coupling of the dual A/" = 4 SYM. Notice that we here 
normalized the R-charge such that a transverse complex scalar $3 in the A/" = 4 has 
+1 i.e. -R($3) = 1. For (l+l)-dimensional fermi liquids, the value of the conductance 
is quantized into ^ times integers (T is the transmission probability), known as the 

On the other hand, in our case, it is quantized in the unit of -y 



Landauer formula 



56 



A 



and we may estimate 7" ~ 4^ for the A/" = 4 SYM media. 

The second solution is a suspended rotating Dl-brane which connects the boundary 
two points P and Q as depicted as the picture (b) in Fig. |2} We describe its proflle by 

V9 = wt + g{r) , 

X = xi = x{r) . (2.47) 

The action looks like 



Sdbi = -Tdi I dtdrVK , (2.48) 
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where 



r^ — iJ^ 



K=' f^+r\r^-u^)x'^ + r^g'^ . (2.49) 

Using the symmetry under constant translations of x and g, we can simply reduce the 
equation of motion into 

9{r) 



2 



x'{r) = , , (2.50) 



where A and B are integration constants. 

Define r^ as a solution to B^{r^ — uj^){A'^r'^ — 1) — 1 = 0. The divergence of x' at r = r^, 
shows that r = r^, is the turning point of the suspended Dl-brane. Thus the Dl-brane 
starts from r = oo and turns around r = r^: and gets back to the boundary r = oo . The 
length L in a; direction of this suspended Dl-brane is 

dr 

(2.51) 



r^^B^ir^ - oo^){A^r^ - 1) - 1 
We can find that g{r) ~ const — ^ + ■ ■ ■ in the r — )■ oo limit. Thus the current fiow is 

^ = % - % = , (2.52) 

dt A A ' ^ ^ 

and the energy fiow is 

^ = ^-^ = 0. (2.53) 

dt A A ^ ^ 



The reason we find that -^ = ^ = is that we employed the simple ansatz (2.47) 
and if we want to violate this we need to consider truly time-dependent solutions which are 
too complicated. The parameters which describe this system are u, A and L. Remember 



that in the previous solution of thermal Dl-branes, we had the constraint (2.9) from the 
regularity of the solutions and the the R-charge current is determined from the total R- 
charge. On the other hand, in the connected Dl-brane solution, they are independent as 
there is no thermal equilibrium even in the fiavor sector. 

Moreover, because r^, > u, all of these suspended rotating Dl-branes do not have any 
horizons in the induced metric 



c?s-„^ = -(r^ - uj^)dt^ + 2ujg'dtdr + ( ^ + r\'^ + g'^ J 



dr^ . (2.54) 



We can compute the time delay along the null geodesic in the induced metric similarly 
to the global AdS case. The induced metric on the suspended solution can be rewritten 
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as follows 

as,„j = — (r — a; j dt , ar + ^„, „ ^, , ,„ „ -r -dr 

*'"' ^ ^\ rV52(r2-u;2)(^2^2_l)_l ^ fi2(^2 _ ^2)(^2^2 _ l) _ l 

= -{r'^ ~uj^)dudv, (2.55) 

where 

« = *-;* ^(-^'■' + -' 



r2^52(r2-a;2)(A2r2-l)-l ' 

v = t+ dr ^ ' 2.56 

' r2^52(r2-cc;2)(y4V2-l)-l 



Thus the time delay of inassless scalar field propagating from a boundary Q to another 
boundary P is given by 

, , B{Ar^-u) r , B{Ar^ + u) 

At = - dr ^ ^ = + dr ^ ' 



oo r2^52(^2_^2)(^2^2_i)_i J^^^ r^ y^ B\r^ - uj^){A^r^ - 1) - 1 

2AB 

dr , . 2.57 

y/B^{r^-uj^){A^r^-l)-l 

Actually, we can show the fluctuations in 6 direction around the Dl-brane are described 
by massless scalar flelds only when A = -. We therefore find 

r dr , , 

At = 2 , (2.58) 

which is much larger than L if i? ^ ^. Thus, as opposed to the global AdS case, we 
can have a large time delay At ^ L and this is due to the strong interactions with the 
deconfined A/" = 4 SYM. 

3 Properties of Thermal Vacua from Rotating D- 
branes in AdS 

In this section we calculate various physical quantities from our rotating D-branes and 
confirm that they are indeed thermalized as we expected from the induced metric. 

3.1 Action of Fluctuations 

First we need to show that the fiuctuations on the rotating D-branes see the horizon which 



we found from the induced metric. Consider a Dp-brane in the background metric (2.16). 
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We choose the gauge where the worldvolume coordinates ^" = x", a = 0...p. The 
transverse coordinates are x^ , J = p + 1 ■ ■ ■ 9. Let us write the metric (3.59) as 



ds^ = gabix", x^)dx°-dx^ + Gij{x"-, x^)dx^dx^ 



Expanding around a classical solution Xq(x°') 



x\x^) = x^(x") + y\x^) 



we find that upto terms of order 0{y'^) the induced metric is 



lab = l^ab + Mab + N, 



ab 



where 7^^ is the induced metric for x^. The matrices M and A^ are given by 

Mab = {digab)y^ + GuidaV^dbx:^ + dby-^daX^) + {dKGij){daxldi,xl)y 
Nab = Gjjday'dby' + {dKGij){day'dt,x:/,+dby'dyo)y'' + 
-{dKdLGij)y^y^daXQdbXo + -{didjgab)y^y 



J„.J 



The DBI action is then 



Sdbi = -TpfdP+'^ Vdet lab 



So + ^ I d^^'^V^o 



(3.59) 



(3.60) 



(3.61) 



(3.62) 



lo'Nab - lifMbcl','M,a + \{l^'Mbay + 0{y') 



(3.63) 



where 5*0 is the action of the background classical solution Xq. 



3.2 Fluctuations around Rotating Strings 

In general the fiuctuations are coupled with each other and have complicated actions. For 
fluctuations around the rotating Dl or Fl solution considered above the action simplifies 



drastically. In the {u,v) coordinates introduced in (2.17), the matrix elements of M and 
N yield the following simple result 



M„ 



0, 



M„ 



2{d,yn{d,^o) 



M„ 



(d^ynid.Vo) 



(3.64) 



and 



A^, 



ab 



Gijday^dby-^ - ida(podb(po){y^ 



(3.65) 
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Plugging this into (3.63) we see that in the first term only the first term of (3.65) survives 
since jQ^'da^Podb^Po = while the contributions from the second and third terms precisely 
cancel each other, leading to the following action for quadratic fluctuations 

Tdi 



S,. 



d\^^oltGij{i\4)day'dby' . 



(3.66) 



In particular, the fluctuations of </?, 9 i.e. all fluctuations in 5''^ directions are minimally 
coupled massless scalars on the worldsheet, while the fluctuations of the boundary gauge 
theory spatial directions x*,i = 1 ■ ■ -3 have an additional factor of r^ coming from the 



fact that Gij = r'^Sij along these directions. 



For a Dl or Fl rotating with a constant angular velocity, it is easy to see that all 
fluctuations perceive a horizon at r = a; and its temperature is Th = ^. For example, 
the action of the fluctuations in x* direction, denoted by X, is given by 



S 



-T. 



Dl 



I 



dtdrJl - :,idrXy + r\r'^ - u'^){X'Y . 



ijj^ 



(3.67) 



3.3 Retarded Green's Functions from Rotating Dl-branes 



As we discussed above, a rotating Dl or Fl corresponds to a source in the (0 + 1)- 
dimensional defect CFT. In this subsection we will discuss the retarded Green's function 
of dual operators to the worldvolume fields. 

For the transverse scalars in the S*^ directions, this can be done using the standard 
bulk-boundary relation |2l|3] in the presence of horizon j8J, since these are minimally 
coupled fields on the worldsheet. Consider for example the worldsheet field if>. This is a 
massless scalar. Let us, however, generalize the analysis and consider a scalar field $ with 



arbitrary mass m which propagates in the induced metric (2.21 ). They correspond to the 



infinitely many stringy excitations of open strings. The equation of motion is given by 



((r^ - Lp\vf)^')' + 2$' - m2$ = . 



(3.68) 



where prime denotes derivative with respect to r and dot denotes derivative with respect 
to V. In particular, we consider the constant angular velocity case, i.e., ^'{v) = uj. Its 
exact solutions with energy Vt look like 



$(w,r) 



^—iVlv 



A 



r + uj 



r 



u 



in 



Pf /'"(r/cj) + B 



r + UJ 



r 



u 



id 



Qf/-{r/u) 



(3.69) 



where k = ^'^+^^ — 1; Pbi^) ^^^ Qti^) ^^^ ^^^ ^^^^ ^^^ second kind of the Legendre 
bi-function. By imposing the ingoing boundary condition at the horizon r = u, we simply 
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find B = 0. The retarded Green's function of the operator O dual to the scalar field can 
be computed from the ratio - when we expand the scalar field as $ ~ ar'^ + (3r~^~'^ in 
the boundary limit r — )■ oo. Finally we thus find that it is given by the ratio of gamma 
functions as follows 

1 [ — K — 1\1/UJ) 

up to some numerical factor independent of Vt. Indeed, for generic values of m, it has non- 
trivial poles only in the lower half plane and has the correct property of retarded Green's 
functions of thermal systems. Particularly for a massless scalar field, the dependence on u 



drops out in (3.70) by chance and the thermal effect is not explicit. However, the same is 
true for the retarded Green's functions in massless free field theory at finite temperature 
and thus we can still conclude that the fiavor sector dual to the our probe D-string is 
thermalized. For example, this issue is clear by comparing the time-ordered two point 



function (A. 142) with the retarded two point function (A. 144) of a harmonic oscillator 
given in the Appendix A. This thermal property will also be shown more clearly in the 
Brownian motion analysis given in the next subsection. 

The fiuctuations in the x* direction have an additional factor of r^ in front of the 
Lagrangian, and leads to the following purely ingoing solution which coincides with the 
one in 



47 



A-(V.)^^(^)'^ (3.71) 

For large values of r, the leading term is a constant, but the coefficient of the sub leading 
1/r term vanishes. In this case, however, the retarded Green's function cannot be simply 
read off from the ratio of the subleading to the leading term, since the equation of motion 
is different, and the analysis of |j8j has to be redone. 

3.4 Brownian Motion of Endpoints for Rotating Dl-branes 

The thermal nature of the state produced by time dependence becomes clear from a 
slightly different calculation - the fluctuation of the end-point of the string. In |47j and |48j 
it has been shown that the fluctuations of a string suspended from the horizon of a AdS 
black brane ended at a flavor D-brane near the boundary of AdS are dual to Brownian 
motion of the corresponding quark in the hot A/" = 4 gauge theory. In this case the bulk 
black brane metric induces a worldsheet metric which has a horizon. The fluctuations 



then reflect Hawking radiation from the worldsheet horizon. The results in 47 and 48 
are consistent with the fluctuation-dissipation theorem. 

22 



In the D-brane solutions considered above, the bulk metric has no horizon. However 
due to the motion of the D-brane, the induced metric on the worldvolume can develop 



a horizon. Since the fluctuations of 47,48 comes purely from properties of the induced 
metric it is natural to expect that a similar phenomenon appears in our case. 

We therefore consider a rotating Dl (or Fl) brane ending on a flavor D-brane situated 
close to the boundary at r = Vc- One example of such a flavor D-brane is the rotating 
D7-brane which will be discussed in the next section and in this case we need to require 
u < Tc- As usual a flnite r^ means that the dual description of the string is a flnite mass 
monopole (or quark), m = TdiTc. Following [47] we now compute the normal ordered 
Green's functions of the end-point of the string, (0, f/| : [y^(t,rc)y^{t',rc)] : |0, f/) in a 
Unruh vacuum state |0, U) of the worldsheet theory. This quantity is equal to the mean 
square fluctuation of the end-point with a trivial divergence due to zero point energy 
subtracted. We will consider the case of a constant spin. 

For constant spin, ifoiv) = u)v the induced metric is 

rfs^ = {-fo)abdCd^'' = "^drdv - (/(r) - uj^)dv^ = -{f{r) - uj^)dudv , (3.72) 

where 

u = v -2r^ , r^= -—- . (3.73) 

J fir) - u^ 

The coordinate r^ is the tortoise coordinate on the worldsheet, where the boundary is at 

r* = and the horizon is at r^ = — oo. The UV cutoff corresponding to the D7-brane 

means that the range of r^ is now — oo < r^ < Ti^^c-, and we need to impose Neumann 

conditions on the fluctuation fleld y'^ at r^ = r^ c. 

The if fluctuations behave as minimally coupled massless scalars in the induced metric 

produced by the rotating string solution. It turns out that for these fluctuations the UV 

cutoff can be removed and we can treat the problem in the full range — oo < r^ < 0. the 

mode expansion is given by 

l/^(t,r,) = / -^^ cos(z/rJKe-^^* + /i.e.] , (3.74) 

J 2Tiy2u 

where 

i^]^{u + v). (3.75) 

Note that on the boundary i = t. The oscillators aj, create "Schwarzschild" particles from 
the Schwarzschild vacuum |0, S) : a,y|0, S) = 0. 

Consider now the case where the background is AdS^ x S^ in the Poincare patch, i.e. 
/(r) = r^. In the Unruh vacuum the correlators of a^y are thermal with a temperature 
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T = 1//3 = uj/2-Ji, i.e. 

(0,f/|4a,,|0,f/) = ^^Ji;^. (3.76) 

Using this it is straightforward to calculate the mean square displacement in the ip direc- 
tion, 

((Ay^(t - t')f) = (0, U\ : [y^it) - y^(t')]' : |0, U) , 

2 f di^sm^{iy{t-t')/2) 



TT J u ef^" -I 

1 /sinh(7r(^-f)/^\ 



Therefore 



i\2 



{{Ay^t-t'))') ~ ^%#, vr(t-t0«/3 



12/3 
((Aynt-t'))^) ~ ^^-^log[2vr(t-t')//3], vr(t-0»/3, (3.78) 

exactly as in Brownian motion. The transition from the short time ballistic behavior to 
the long time diffusive behavior happens at a time scale set by the inverse temperature, 
which is the only scale in the problem. The result is identical for the y^ fluctuations, as 
well as fluctuations in the S^ directions. 

The fluctuations in the x directions are not conformally coupled scalars. The equation 
is in fact identical to the equation for fluctuations of the end point of a string in the 



background of a BTZ black hole, as in 47 . In this case, the UV cutoff at r = Tc is 
essential. The final result is |47j 

{[Ay\t-t')r) ~ ^^^, t«m/3^ 

~ I3\t-t'\ , t>m/32 . (3.79) 

Unlike fluctuations in the internal 50(6) space, Brownian motion in the physical space on 
which the gauge theory lives is present only if the mass of the monopole or quark is finite. 
This is of course expected. Furthermore the time scale which controls the transition from 
a ballistic behavior is not simply the inverse temperature, but involves the mass in an 
essential way. 

As argued above a rotating Dl-brane corresponds to a time-dependent coupling in the 
defect (0 + l)-dimensional CFT. Brownian motion is a manifestation of thermalization of 
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this defect CFT, or equivalently the hypermuhiplet sector of the theory. This Brownian 
motion is a physical consequence of the Hawking radiation on the world-sheet and this 
corresponds to the thermalization of the hypermultiplet sector (or defect sector). Notice 
that this is not directly related to the Hawking radiation in the bulk AdS which conies 
from the thermalization of the A/" = 4 super Yang-Mills sector. Indeed their temperatures 
can be chosen differently. 

3.5 Conductivity from Rotating D3-branes 



For rotating Dp-branes of section (2.3), there are several other response functions which 
carry signatures of thermalization. A good example is a probe D3-brane configuration ob- 
tained by taking T-dual of the rotating Dl in x, y direction. This is a non-supersymmetric 



configuration in AdS^, x S^. The rotating D3-brane solution is obtained from (2.24) and 



(2.25) by setting p = 3. Its induced metric is given by the black hole geometry (2.26) at 



p = 3 and its Hawking temperature is Th = ^(^■ 

The worldvolume can now contain a nontrivial gauge field, and the DBI action of the 
D3-brane looks like 

C = Jr^ + rV^ - ^202 - (1 + r^ip'^)Fl - r^Fl + (r2 - ip^y^F^ + 2r'^Fr^Ft^ipip' , 



^6 



:F4 - r^Fl + r\r^ - uj^)Fl , (3.80) 



This shows that the gauge field is propagating in a space-time with a horizon. 

We would like to calculate the electric conductivity in the probe D3-brane background 



by applying the method developed in AdS black holes [8, 11 . We regard the abehan 
gauge field on the probe D3-brane as an external gauge field for which we define the AC 
conductivity cr(z/), where z/ is the frequency. In other words, here we are thinking the 
charge with respect to the global flavor symmetry as the electric charge to define the 
conductivity. 



The DBI action (3.80) can be expressed using the induced metric Qind- 



F^,F<-P 



Ld3 = V-det(7,w\/l + ^-^ . (3.81) 

It follows that the equations of motion are 

pal3 

, 1=0. (3.82) 

^1 + F2/2 ' 
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Before computing the conductivity, we have to find a configuration of gauge field with 
non zero gauge potential Aa = {At{r) = $(r), 0, 0, 0) which obeys the following equation: 



dr 



„2 . / r'^+r^uj^+ui'^ ^/ 



'1 r^+r^i^^+cj-^ ^/2 



(3.83) 



We can easily find a solution with the boundary condition $(r) = at r = cu: 

r,2 



^(r) 



ds- 



The gauge potential behaves near r ~ oo like 

$(r)=^-^ + ... , 



(3.84) 



(3.85) 



and (/i, p) stand for the chemical potential and the electric density, respectively. Accord- 



ingly, we can compute the electric density using (3.84) as follows: 



p = lim <!>'{ry = C-^/2 _ 



(3.86) 



Now we consider a small perturbation A^ = Ax{r)e "^^^ on this background. We use 
z = 1/r coordinate for simple calculation in the following. The equation we solve is that 



^ . y, _ z\Qz^ + z^ + z' -2) + Cz{z^ + z' + A) ^^ 



y'A^ 



{l-z^){l-z^) 



, (3.87) 



where we rescaled z — )■ z/u, and define C = Cuj^ and z/ = p/u. The gauge field is an 
ingoing wave near the horizon z = 1, and it thus can be expanded as follows: 



A,{z) = (1 - z)-'^v-3(l + ai(l -z) + a2(l - zf + . . . 
while it behaves at the boundary z ~ as 

A^{z) = A^^^+A(^h + ... . 
We can obtain the conductivity 



(3. 



(3.89) 



a{u) 



iA 



(1) 



uA'i 



(0) 



(3.90) 



In Fig. Hwe plot the conductivity obtained by solving (3.87) numerically. The behavior 



that Recr(z/) approaches to a finite constant is known to common in (2+l)-dimensional 
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critical theories 11 . As opposed to the purely bulk calculation of the conductivity from 



the AdS charged black holes, where we have Recr(z/) oc 5(z^), instead we observe a smooth 
Drude-like peak both in Recr(z/) and Imo"(z/) near z/ = in Fig. |4J This is because our 
probe approximation introduces dissipation into the bulk and has the finite DC conductiv- 



ity. A very similar behavior has been already obtained in 29 , 33 from the probe analysis 



of D-branes in AdS black holes. Here we managed to get a similar physics by using a 
rotating D-branes in pure AdS. 



Re(cr) 





Figure 4: The frequency dependence of the conductivity. Four lines are depicted for 
p = 10, vio, 1 and 1/vIO, respectively when to = 1. The real part of the conductivity 
approaches to Recr(i^) = 1 line as p decreases, while the imaginary part becomes zero. 



3.6 Comments on Entropy for Rotating D-branes 

Before we finish this section, we would like to comment on a possible definition of entropy 
of rotating D-branes. Since we find a thermal horizon in the induced metric, it is natural 
to expect that there exists non-trivial entropy in such a system in spite of the absence of 
real black holes. 

There is one natural candidate of such an entropy for probe Dp-branes: it is calculated 



from their free energy by Wick rotating the time coordinate r in (2.26). Consider first 
a rotating Dl-brane in AdS^ x S^. The free energy is simply calculated from the DBI 
action as follows 

poo 

F, = TdiR^ / dr , (3.91) 

where we revive the AdS radius R for convenience. The entropy is obtained by taking 
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derivative w.r.t. Th 



2lT 



dFi ^_^2. 




^1 = - ^ = 27rR'Tni = 27r J -- = 2 v^^ . (3.92) 

The linear dependence of the rank A^ of gauge group agrees with what we expect for the 
flavor contribution in the planar limit. 

Let us move on to the higher-dimensional D-branes. As one such example, we would 
like to consider the rotating D3-brane in AdS^ x S^ discussed in the previous subsection. 
Then its free energy looks like 

Fs = T^si^Vs / dr 



4 



= TD3R''V2('f-^r^co' + ku'\ , (3.93) 

where k ~ 0.135 is a numerical constant; V2 is the area of non-compact space defined by 
the intersection of the D3 with the AdS boundary. We also kept the UV cut off roo(— )■ 00) 
explicitly. By taking the derivative with respect to the temperature Th = ^^^, we get 
the following entropy 

OF, 27r 4 2, 2x 2NV2 [Th r^A .^ Q,^ 

= ~Qf~ = 7^^^3^ L [r^uj -3buj ) = -^— ( — - kTjjj , (3.94) 

where a = 1/roo is the lattice spacing (or UV cut off) of the dual CFT and k = 27r-\/3A; ~ 
1.47 is a numerical constant. Thus this entropy is UV divergent and this behavior is 
always true for rotating Dp-branes with p > 1. One may think this is unnatural because 
the entropy which we calculated in this way should be the thermal entropy of the hyper- 
multiplet, which should be finite. This might suggest that the above way of computation 
of entropy via our naive Wick rotation is not correct. However, we would like to propose a 
possible interpretation of this divergent entropy. Recall that in our probe calculation, the 
interactions with the A/" = 4 SYM sector were already incorporated as the supergravity 
background. This means that the A/" = 4 SYM sector was already traced out in the den- 
sity matrix. Therefore the entropy which we calculated in the above should be regarded 
as the entanglement entropy between the flavor sector and the A/" = 4 SYM sector. In 
thermal equilibrium, we can neglect the entanglement. However, in our non-equilibrium 
systems where the two different temperatures coexist, there may be non-trivial entangle- 
ment entropy which can be UV divergent. Notice that the UV divergence of the standard 
entanglement entropy is holographically explained by the infinite volume of AdS space [57] 
and the origin of the divergence here looks similar. It is a quite intriguing future problem 
to confirm this from field theory calculations. 
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4 Rotating D7-branes in AdS^ x S' 



So far we have assumed that the rotating Dp-brane is point-hke in S'^. However, this 
class of D-branes misses important supersymmetric configurations because the latter often 
wrap on S"' {n > 1) in S'^. One such example is the probe D7-branes in AdS^ x 5*^, which is 



dual to the A/" = 2 flavor hypermultiplets coupled to the M = 4 super Yang-Mills 26 - 28 
Notice that in this case the space-time of the dual field theory is (3 + l)-dimensional even 
for the probe D7-brane sector. Motivated by this, in this section we will study the rotating 
D7-branes in AdS^^ x S^. This system has two independent parameters i.e. the distance 
m between the D7-brane and the background D3-branes and the angular velocity u. This 



setup has been already studied in 49 , 50 . Though our D7-brane solution presented below 
is essentially the same as the one in |J4^, we will study the solution from a different 
viewpoint focusing on the presence of the thermal horizon in the induced metric. 

4.1 D7-brane World- Volume 

The D7-brane world-volume coordinates are extending in the AdS^ and S^ included in S^ 
with a rotation in ip direction. This is specified by 

9 = 9{r), (p = ujt + g{r). (4.95) 

The DBI action now becomes 

C = -r^ cos^ eVZ , (4.96) 

where 

L = l + r^ sin^ 9^'^ - ^' ""'f ^ + r^9'^ _ uj^ sin^ 99" . (4.97) 



The equation of motion for (f reads 

d f r^ cos^ 9 sin" 9ip' 



dr \ y/Z 

The equation of motion for 9 is given by 



. (4.98) 



2^/2 ui^ , ,2/1/2 



'ir' cos^ 9 sin 9^ - r^ cos^ 9 sin ^^^ ^ + dr 



\ -^ A^"^ — w^sin^'^''^' 
r^cos'^6'^ 



Vl 



= . 

(4.99) 



The solutions to the equation (4.98) should be either 

g\r) = , (4.100) 
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or 



g'ir) 



r^ sin 6 
where A is an integration constant. 







(4.101) 



Let us first study the case (4.100). When w = 0, we can find the solution rsin^ 



771. 



This corresponds to the D7-brane separated from D3-branes by the distance m, which 



is dual to the mass of hypermultiplets. Now consider solving the equation (4.99) with a 
boundary condition 6 = 7i/2 at r = a where a is an arbitrary constant. The function 
r{6) for 6 = i\ j1 + b for an infinitesimally small b is then determined by the equation of 



motion (4.99) 



r(7r/2 — (5) ~ a + 



4a^ - 3cj^a 



52 + . 



(4.102) 



Using this we can solve for the function riff) numerically for < w < a. In this case the 
solutions does not differ qualitatively from the a; = solution and there is no horizon, as 
can been from the left graphs in the Fig. |5| If we increase lo such that u; > a, we cannot 



obtain smooth solutions from this ansatz (4.100) 



r sin(«) 
2.0 - 




/ / 



I ^.*— . 



0.0 0.2 0.4 



0.0 0.2 0.4 0.6 0.8 



i.O 1.2 1.4 



Figure 5: The shapes of the rotating D7-brane with w = 1. In the left panel where we 
consider f?' = case, we depict the curves for a = 1.01, 1.1, 1.2 and 1.5. In the right panel 
where we consider g' 7^ 0, we depict them for a = 7r/2, 7r/2.1, 7r/2.5, 7r/3, 7r/4, tt/IO. The 
dotted regions are inside the horizon. 



Next we would like to consider the case when (4.101) is satisfied. The denominator 
A'^r^ cos^ 9 sin^ ^—1 can become negative. To avoid the resulting singularity in the function 
L, we need to require that r^ — u^ sin^ 9 and A'^r^ cos^ 9 sin^ 9 — 1 vanish at the same time 
i.e. 



a 



uj^ sin^ a 



Al'o^ cos^ a sin^ a 



(4.103) 
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where we set r = a and 6 = a at this point. This point will turn out to be the horizon. 



(4.104) 



If we expand the solution near this point, we find from the equation of motion (|4.99|) 

6(r) ~ a + /3(r — a) + ■ ■ ■ , 
where the (negative) constant /3(< 0) is determined by 

3 + 2Pu cos a - 3(3^uj^ sin^ a = . (4. 105) 

The solution to this equation for /3 is 



/3 



(cos a — V cos^ a + 9 sin^ a) 



(4.106) 



3uj sin a 
Thus if we input the values (a, u) we can uniquely find the solution numerically by solving 



(4.99) starting from the thermal horizon point r = a. They are shown in the right graphs 



of Fig. pi Notice that the two solutions (4.100) and (4.101) of the D7-brane match at 
u = a (or a = vr/2) as depicted in Fig. 16} 




Figure 6: The D7-brane is shaped like almost fiat brane for u < a, and becomes singular 
when u = a. Then it always have a horizon on its worldvolume for u > a. 



4.2 Induced Metric and Thermal Horizon 

The induced metric reads 



ds' 



(r^ - w^ sin^ e)dt'^ + 2u sin^ eg'{r)dtdr + cos^ Odnj 
+ {\ + e'^ + sin^ 9g'A dr^ + r\dx^ + dy^ + dz^] 



After the coordinate change t = t + h{r) with 

uj sin^ 6g' 



h\r) 



r^ — u 



2sin2^ ' 



(4.107) 



(4.108) 
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we can rewrite the metric as follows 



\r^ r^ - uj^ sm 9 

+ r'^{dx^ + dy"^ + dz^) + cos^ OdVtl . 



1 



r^ sin^ 6g''^ 



(4.109) 



Thus the horizon is identified with r = wsin^ as already mentioned. The topology of 
horizon is given by R^ x S^. 

The Hawking temperature can be found from this metric 



T 



u /sinQ;tana;(4sina — /3a;(3 + /SwcosSo;)) 



27r 



1 + (3'^uj'^ sin a 



(4.110) 



We can see that T is monotonically increasing function of a when we fix cj. At a = w, we 
find T = oo, while at a = 0, we have T = 0. The fiavor mass m = linir^oo f sin 6' of this 
solutions is plotted as a function of a in the left panel of Fig. [7} The dependence of the 
temperature on the mass is also presented in the right panel of Fig. [7j 





Figure 7: [Left] The dependence of the function m = limj._i.oo ^ sin ^ on a in the second 



solution (4.101). We set a; = 1 and connected the two solutions (4.100) and (4.101) at 



a = 1. [Right] The dependence of the temperature on the mass m at a; = 1. For m > 1.15, 
there is no horizon. 



4.3 Dual CFT Interpretation 

When the D7-brane does not rotate, our setup is the same as the standard D3-D7 system. 
The dual CFT is described by the four-dimensional A/" = 2 SU{N) SQCD which is defined 
by the A/" = 4 SYM coupled to a A/" = 2 fiavor hypermultiplet in the large N limit. The 



time-dependence oi (p{t) = ut (4.95) at the boundary r = cxd is dual to the time- dependent 
mass 



m{t) 



me 



iip{t) 



me 



iuit 



(4.111) 
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In terms of the superpotential we have 

W = Tr[$i, $2]$3 + Q^sQ + m{t)QQ , (4.112) 

where $i,$2,$3 are the transverse complex scalar fields on the A^ D3-branes. {Q,Q) 
are the hypermultiplets from D3-D7 open strings. We will also denote the fermions in 
the hypermultiplet by {TpQ^ipq). The original SU{4) = S0{6) R-charge of A/" = 4 super 
Yang-Mills theory is now broken to 

^0(4) X U{1)r = SU{2)r X SU{2)l x U{1)r . (4.113) 

The S'0(4) rotates the R^ described by $i and $2- U{1)r corresponds to the rotation 
of the remaining R^ i.e. $3. The R-symmetry of A/" = 2 super- Yang Mills is given by 
U{2)r = SU{2)r X U{1)r, If we denote the spins under SU{2)r x SU{2)l x U{1)r by 
{Jr,Jl, ^[/(i)); the complex scalar fields {Q, Q) in the hypermultiplet belongs to (1/2, 0, 0). 
On the other hand, fermions {iPq,iPq) belong to (0,0, ±1). The transverse scalar fields 
$1, $2 belong to (1/2, 1/2, 0), while $3 to (0, 0, 1). 



Now we would like to come back to the superpotential (4.112). After integrating the 



fermionic coordinates of the superfield, we find the time-dependent bosonic potentia]^ 



g(<l>3 - me^'^*)(<l'3 - me''^^)Q + Q($3 - me^'^*)(<l>3 - me'^^) Q . (4.114) 
as well as the time-dependent fermion mass term 

dx^ me'^'ipQipQ + (h.c.) . (4.115) 



The AdS/CFT claims that the above time-dependent system will lead to a thermal vac- 
uum. These time-dependent potentials triggers a physics similar to that in quantum 
quench and strong interactions in the hypermultiplets will be expected to lead to a ther- 
malization. However, since the hypermultiplet is coupled to the A/" = 4 SYM sector, the 
energy eventually will dissipate, as we have explained in the simpler model of rotating 
Dl-brane previously, schematically summarized in Fig. [l] It will be a quite interesting 
future problem to confirm this prediction from the field theory calculations. 

^ One may notice that if we redefine $3 as $3 = $36*'^*^*-', then the system is equivalent to the standard 
flavored N =2 SYM with the R-charge chemical potential At — lu only for the field $3 without any time- 
dependence. However, from this viewpoint we need to take a time-dependent state and its time evolution 
will be described by this time-independent Hamiltonian. This is because we know from AdS/CFT that 
the hypermultiplet is initially exited and has non-zero R-charge instead of the M — A SYM sector. 
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5 Quantum Quench and Thermalization via AdS/CFT 

So far we have mainly studied stationary configurations of the Dp-branes. We would 
now like to discuss truly time-dependent D-branes in AdS spaces and its interpretation via 
AdS/CFT correspondence. The goal of this section is to holographically model quantum 
quench - the time evolution of a system following a sudden change of a parameter in 
a given quantum field theory, such as the mass of scalar fields. The phenomenon of 



quantum quench (see e.g. 16,17,19 22 for approach from conformal field theories) has 
been intensively studied in condensed matter theory recently. 



5.1 Quantum Quench in Free Field Theory 

Some interesting aspects of the effect of quantum quench can be in fact observed in a free 



quantum field theory 19,21 . Let us consider a d-dimensional massive bosonic free field 



theory with Hamiltonian 



H 



(Tk 



-VTfcTr^fe + -Uf,(f>k(p~k 



(5.116) 



with a dispersion relation ul = k'^ + m?. 

Each momentum mode is independent of each other, and the propagator is just that 
of a single harmonic oscillator. We therefore consider a single harmonic oscillator with 
momentum k that initially lies in a thermal state with temperature /3o, and then we 
quench its frequency from wofc to Uk by changing mass from itiq to tti at t = 0. At a later 
time t after the quench, the scalar field is given by (pkif) = 0fe(O) coscokt + TCkiO) sinukt/cok- 



The time ordered correlator will be 21 



Cp,ik;ti,t2) =T{Mti)(l>^kit2)) 



2uk 



g-iajfc|ti-t2| _|_ 



^Ok 



1 1 



w, 



Ofe 



uJt 



coth 



Po^Ok 



2ujk 



+ 



UJ, 



Ok 



UJt 



coth cos Uk (ti + ^2) 



cosa;fc(ti - t2) 

(5.117) 



This result (5.117) directly follows from the two point function (A. 146) in a harmonic 



oscillator as explained in the Appendix A of the present paper. 



In real space, the correlator is obtained by taking the Fourier transform of (5.117) 



The important fact is that the third term in (5.117), which violates time translation 



invariance, turns out to decrease and eventually vanish after a long time under quite 



general conditions 19,21 . The time one has to wait to see this effect depends on the 
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wave-number k and increases as \k\ decreases. Notice that the density matrix of this time- 
dependent system is always pure since we start from a pure state at zero temperature and 
it evolves in a unitary way under the time- dependent Hamiltonian. Thus to find a thermal 
behavior we need a process of coarse graining. 

In this way, we can ignore this part and compare the others to the thermal propagator 



with temperature /3 (see (A. 142) in Appendix A for an elementary derivation) 



G,{t,t,M) = -Le— l..-'.l + '°7^-'^' . (5.118) 

2uk Wfc^eP^fe - 1) 

We can read off the effective temperature by equating the coefficient of cosa;fc(ti — t2) 

It follows from this result that the effective temperature depends on each momentum 
mode. 

It is clear that an effective temperature is obtained even when the initial state is the 
vacuum state rather than a thermal state. In this case 

/3.,,,o(A:) = -log/^^^±^. (5.120) 

(^k \^k — ^Ok\ 

It is easy to extend the above calculation to a situation where the mass starts from 
a value mo at t = 0, changes to a value m ^ tuq for < t < T and then finally changes 
back to m = ttiq for t > T. Interestingly the late time correlation functions for t ^ T are 
again thermal with a fc-dependent effective temperature. This is consistent with the fact 
that once the system thermalizes, it remains thermal. 

5.2 Exact Toy Model from Dl-branes 

The framework considered in this paper can be used to construct a very simple toy model 
for a quantum quench in a strongly interacting theory. This is provided by the exact 
time-dependent solution of Dl-brane whose induced metric is given by the AdS Vaidya 



metric studied generally in (2.21). In our current example, its two-dimensional induced 
metric looks like 

ds'^ = 2drdv - (r^ - yD'{vf)dv^. (5.121) 

To describe a quantum quench we can for example take {v = t — 1/r) 

cp(v) = ipo{l + tanh kv) , (5.122) 
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which approaches a step function in the large k hmit. Then the apparent horizon is 

situated at 

r = ^\v) = ^^ . (5.123) 

cosh kv 



In this model, the thermalization is occurred due to the time-dependent coupling (2.15) 
as we discussed before. 

When the apparent horizon is slowly evolving, it is possible to define a local tem- 
perature by considering the analytic continuation of outgoing modes across the apparent 



horizon 58 . In our case this leads to a "local" temperature given by ^'{v). While this ex- 
pression is not strictly vahd when the evolution is fast, we can still use this to get a rough 
estimate of the temperature. In an eikonal approximation, the Hawking particles emitted 
at some point (r, t) on the apparent horizon will travel to the boundary along a line of 
constant u = t+ 1/r and will therefore reach the boundary at time t = u — 1/r = v + 2/r. 
Therefore the temperature perceived at the boundary at time t is given by 

TH{t) = ip\v), t = v + ^ (5.124) 



where we need to first express v{t) using the second equation of (5.124) and then substitute 
this in the first equation to obtain TH{t). 



For a ^(v) of the form (5.122) this is plotted in Fig. M The temperature is highest 
at t = v^'(O) = ipo and decreases to zero at t = oo. Notice that the late time behavior is 
given by 

Tnit) ~ ^ . (5.125) 

This is a consequence of the conformal invariance of our D3-D1 system. 

Quantum quench calculations in the literature have often been done in a closed envi- 
ronment. In such a case, a system got excited suddenly at some time by an external force 
and will reach an equilibrium which are static. On the other hand, in our system we find 
that the effective temperature decreases to zero finally. This is because the hypermultiplet 
sector, which are initially excited, is strongly coupled to the A/" = 4 SYM sector and thus 
the energy will dissipate. 

5.3 D3-D5 System with Time-dependent Mass 

Here we would like to turn to a more interesting example: a holographic dual of quantum 
quench induced by time-dependent masses of scalar and spinor fields. For this purpose 
we consider a probe D5-brane wrapped on AdS^ x S^ in AdS^ x S^. This setup is super- 



symmetric and was first considered in 25 . It is dual to the four-dimensional J\f = A SYM 
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Figure 8: The time-dependence of the effective temperature Tnit) in the Dl-brane model 
under an eikonal approximation. We set k = 4 and (po = 1. 



coupled to a three-dimensional defect CFT. If we write the metric of AdS^ x S^ as 

3 '2 



dr 



ds^ = r\-dt^ + Y^ dxi) + ^ + de^ + cos^ Od^l + sin^ OdQ. 



2 
2 ) 



(5.126) 



i=l 



then the D5 is wrapped on VL2 and (t,Xi,X2,r) directions, while it is trivially localized in 
VL2 and X3 direction. However, in actual numerical calculations, it is useful to perform the 
following coordinate transformation 

1 



T] = r sin 9 , z 



r cos 6 



(5.127) 



for which the AdS^ x S^ metric looks like 

3 7 2 



ds' = iv' + z-')i~dt' + J2dx^) + ^^^ + 



dr]'' , dz' , dnl , rj'^dn'l 



+ 



+ 



i=l 



rf + z~^ z^ + z^if' 1 + z^rf rj^ + z 



(5.128) 



In this new coordinate, we can specify the profile of the D5-brane by the function rj 
■r][t,z). Its DBI action reads 



S = -TD5V4 / dtdz-Jl + zV' - 



^2 



z'* \l ■ (r^^ + z-2)2 

where V4 is the (infinite) volume of the 5*^ x R'^. Its equation of motion looks like 
d f rj' \ Of fj \ 2iri'rf 



(5.129) 



dzy^/Kj dt \(i + z^rj^y^/K J z\z~^ + r]^)3^/K 
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0, 



(5.130) 



where K = l + z^t]'^ - (^^2^^,-2)2 ■ 

When Tj can be treated as a small perturbation so that it satisfies 

z^T]' <C 1 , and z'^rj <^ 1 , 



(5.131) 



we can approximate the equation of motion (5.130) by the simple wave equation i.e. 
7]" — 7] = 0. Therefore, if we require that the defect sector in the dual CFT has the 
time-dependent mass m{t), we can find the simple solution under this approximation 

r]{t, z) = m{t ~ z) . (5.132) 



However, in the interested regions where horizons appear, the condition (5.131) is in- 



evitably violated. This is clear from the expression of Qu component in the induced 
metric 



9tt\ind = -{V + Z ) + 



T] 



(5.133) 



^2 _|_ ^-2 

which should change its sign at apparent horizons. Since in our particular setup, the area 
of a codimension two surface for fixed values of t and z is always only a function of z 
(proportional to z'"^), the apparent horizon for our induced metric, which is defined by 
the vanishing area change along null geodesies, actually coincides with the points where 
gtt vanishes. 



Thus we numerically integrate the equation of motion (5.130) under the two boundary 



conditions presented below. The first boundary condition is imposed at the AdS boundary 
with a UV cut oS z = zq <t^ 1: 



r]{t,Zo) =m{t- Zq) , 
The second one is at a initial time 1 = 1^. 

ri{to,z) = m{to- z) , 



dri{t, z) 



dz 



dr]{t, z) 



z=zo 



dt 



t=to 



-m{t - Zq) . 



m{to - z) . 



(5.134) 



(5.135) 



Below we would like to analyze a specific choice of the mass function m{t) defined by 

m{t) = mo(l + tanh H) , (5.136) 

which describes a sudden quench for large k. 



In both of the above conditions we employed the approximated solution (5.132). For 



the first one (5.134) this is clearly justified because zq is infinitesimally small. On the 

1 



particular case (5.136) 



other hand, for the next one (5.135), the condition (5.131) requires to ^ 



/mok 



m our 
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Figure 9: The function rj{t,z) for the choice k = 10, mo = 0.1, to = 0.01 and Zq 
We plotted the region defined by 0.01 < t < 2 and 0.001 < ^ < 3. 



0.001. 



Following this strategy, we plotted the function ri{t, z) in the Fig. ^for the particular 
choice of the parameters k = 10,mQ = 0.1, to = 0.01 and Zq = 0.001. As is obvious from 



this graph, the kink, generated by the sudden mass change (5.136) in the UV, propagates 



at the speed of light into the IR region. This is qualitatively similar to our previous exact 
Dl-brane solutions discussed in the previous subsection. To show the existence of the 
apparent horizon, we need to examine the sign of Qu of the induced metric. We plotted 
—gtt in Fig. [lOjonly when it is positive. We can confirm a hairpin-shape boundary near the 
line t = z, where gu vanishes. Therefore we can conclude that our time-dependent mass 
at the boundary generates a thermal horizon after a certain time (estimates as 



/mok ' 



Its temperature from an observer at the boundary decreases as time passes and eventually 
goes to zero due to the dissipation of energy into the bulk. From the dual CFT viewpoint, 
this clearly describes the process of the thermalization in the open system when we change 
the mass of scalar fields and fermions at some time. 

Finally it is instructive to compare our result with that of free field theories. In the free 
scalar field theory, after a sudden change of scalar mass, the system should reach a sort 
of thermal equilibrium state. However, as we have reviewed briefly, different momentum 



modes have different temperatures (5.119) 19,21 . This momentum dependent tempera- 



ture arises because the different momentum modes are not mixed with each other in the 



free field theory. A similar issue has been known to occur in integrable systems 18,21 



where the final equilibrium state is described by the generalized Gibbs ensemble, due to 
the infinite number of conserved quantities. However, our system is a strongly coupled 
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Figure 10: The plot of —gtt{t, z) of the induced metric for the same choice as the previous 
plot. We concentrated on the region defined by 0.5 < t < 1.5 and 0.8 < z < 2. We 
plotted the graph only when —gu > and thus the hairpin-shape boundary in the middle 
represents the positions of the apparent horizons. 

quantum field theory as we work in the supergravity regime and we expect substantial 
mixings between states with different momenta. Indeed, our holographic analysis shows 
that the effective temperatures for all momentum modes in the dual CFT are originated 
from the single apparent horizon at a fixed time. Therefore our holographic description 
predicts that the final temperatures for each momenta at the equilibrium should be the 
same in strongly coupled field theories. 
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A Thermal Correlation Functions and Mass Quench 
in Free Theories 

A.l Thermal Correlation Functions 

Here we would like to give an elementary derivation of thermal correlation functions for a 
harmonic oscillator and a free scalar field theory. We start with the standard Hamiltonian 
of a harmonic oscillator 

H=^p' + ^x' = coha+^^ , (A.137) 



where 



t ^ I r-.- ■ P \ _ 1 / /"T- , • P 



a' = —pz I \/uJx — i—p= , a = —= ( JZjx + i^i= . (A. 138) 

v^V v^/ V^V v^/ 



The vacuum |0) is define by a|0) = and we can define the normalized number state 

\n) = -7=f(a''^)"|0) as usual. 

tl 



Using the commutation relation [a, a'] = 1, we obtain 



x(t) = -L (ae-*^* + a+e^^*) , (A. 139) 



The finite temperature density matrix for the inverse temperature /3 is give by 



2 sinh ^ 



P = e-^- = ^^^ ^^, |^o-;r^v^l- (A-140) 



We can calculate various two point functions as follows. First of all, we obtain the 
Wightman function 

posh (^ - iboA 

(mm) - Trie^'^mm] = 7 ^ ' , ,T' ■ (a.mi) 

Zoj smn -^ 

Next, the time ordered two point function is given by 

^, / N / Nv cosh ( !^ — iu\t\) 1 . ,,| coswt 

T{x(t)xm = ^^ ^^ = — e-*"^!*! + —— . A.142 

^ ^ ^ ^ '^ 2a;sinhf^ 2co co{eP-^ - I) ^ ' 

It is also useful to calculate the commutator 

([x(t),x(0)]) =Tr[e-^-^[£(t),x(0)]] = --sinwt , (A.143) 
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which is actually independent of the temperature. This leads to the retarded Green 
function 



Gr = {[x{t),m])R = e{t)Tr[e-^''[x{t), 



X 



u 



-6(t) smut 



(A. 144) 



In the free scalar field theory (5.116) with a mass m, if we work in the momentum basis 



k, the calculations are identical to Harmonic oscillators with the fc-dependent frequency 
Uk = Vm'^ + k^- Therefore the time ordered two point function is found to be (5.118) as 



is clear from (A. 142). 



A. 2 Mass Quench 

Now we would like to analyze the setup of the quantum quench due to the sudden change of 



frequency a; in a harmonic oscillator (A. 137) as considered in 21 . Let us assume for t < 



the frequency is given by ojq and the system is in a thermal equilibrium at temperature 
To = l//3o- At t = 0, the frequency suddenly is changed into w by a certain external force 
and we are interested in physics when t > 0. Since the x and p are continuous at t = 0, 
the creation and annihilation operators change at t = from (ao, %) to (a, a^) via the 
following Bogoliubov transformation 



cosh C cto + sinh C a^ 



sinh C cto + cosh ( Oq 



cosh^ = 



OJ 




+ 



Uo 




(A. 145) 



Using (A. 145) and (A. 139), we can obtain the time-ordered two point function after the 



quench ^1 



T{x{U)x{t2)) = T [Tr[e-''«^°x(ti)a;(t2)]] 



+ 



1 

Wo 



g-«^l*l-*2| _|_ 



CJo 

4 



u^ 



UJr 



coth 



/3oWo 



1 

%0 



cosuiti 



OJa 



U' 



coth cosa;(ti + ^2) 



(A. 146) 



In section 5.1 we discusses the quantum quench in the free massive scalar field theory, 
where the mass is suddenly changed from ttiq to m. In this case we can simply replace 
u and Uq with the fc-dependent frequency Uk = y/W^rm^ and Uq^ = \/k'^ + mg, respec- 



tively. This leads to the two point function (5.117). 
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